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Cosmic microwave background data shows the observable universe to be nearly fiat, but 
leaves open the question of whether it is simply or multiply connected. Several authors 
have investigated whether the topology of a multiply connect hyperbolic universe would 
be detectable when 0.9 <Q < 1. However, the possibility of detecting a given topology 
varies depending on the location of the observer within the space. Recent studies have 
assumed the observer sits at a favorable location. The present paper extends that work 
to consider observers at all points in the space, and (for given values of and Oa 
and a given topology) computes the probability that a randomly placed observer could 
detect the topology. The computations show that when f! = 0.98 a randomly placed 
observer has a reasonable chance (~ 50%) of detecting a hyperbolic topology, but when 
n = 0.99 the chances are low (< 10%) and decrease still further as C approaches one. 
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1. Introduction 

Analysis of recent cosmic microwave background (CMB) data suggests an approxi- 
mately flat universe with the^total energy density parameter almost surely lying 
in the range 0.9 < il < I.I.EId The near flatness of the observable universe does not 
preclude a multiconnected spatial topology, but may push the topology to a scale 
larger than the horizon radius, making it difficult or impossible to detect. Recent 
studies have examined the extent to which a nontrivial topology may oxjnay not 
be observable in a locally spherical universa3 or a locally hyperbolic oneJj'BQO In a 
locally flat universe there is no a priori relationship between the topology scale and 
the horizon scale, so a great deal of luck would be required for the two to coincide. 

In their most recent study of multiconnected hyperbolic universesjj Gomero, 
Rebougas, and Tavakol examine the seven smallest known hyperbolic topologies 
and find that for a set of cosmological parameters given by Bond et al% with 
O = 0.99, five of the seven topologies would be potentially detectable using CMB 
methods, while for a set of parameters given by Jaffe et aZ.i with = 0.98, all 
seven would be potentially detectable. 

A topology is considered potentially detectable by an observer at a point p in the 
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space if the observer's horizon radiu^ rhor exceeds the injectivity radius^ rinj (p) 
at p. In a 3-torus and in many spherical topologies, the injectivity radius Vinj is 
constant throughout the whole space, so if the topology is potentially detectable 
by an observer at point p, it is detectable by any other observer at any other point 
q in the same space. In hyperbolic topologies, by contrast, the injectivity radius 
varies from point to point, so a hyperbolic topology might be detectable by an 
observer at point p (where the injectivity radius ri„j{p) is small), but undetectable 
by a different observer at some other point q (where the injectivity radius ri„j{q) 
is large). 

In their most recent work,Ll Gomero, Rebougas, and Tavakol consider the ques- 
tion of detectability from the most favorable point in the space, that is, from a point 
of minimal injectivity radius. The present article extends their work by considering 
the detectability of the topology at arbitrary locations in the space. The variation 
of the injectivity radius across the space will be summarized in an injectivity profile 
showing what fraction of the space's volume has a given injectivity radius. Com- 
bining the injectivity profile (Section ||) with an estimate for the horizon radius 
(Section ||) reveals the probability that a randomly placed observer could poten- 
tially detect the topology, that is, it tells the fraction of the manifold's volume in 
which rhor > finj (Section 



2. Injectivity Profiles 
2.1. Definition 

An injectivity profile is a histogram showing how much of a manifold's volume 
has a given injectivity radius. For example, in the simple histogram in Fig. 1 the 
first bar shows that the injectivity radius lies in the range [0.2, 0.3] for 10% of 
the manifold's volume; the second bar shows that the injectivity radius lies in the 
range [0.3, 0.4] for 20% of the manifold's volume; and so on. In the studies presented 
below, each histogram will have 240 bins of width 0.0025, spanning the range of 
injectivity radii from to 0.6. In the limit, as the bin width goes to zero, the 
assignment of a finite volume percentage AV/V to each bin of finite width Ar 
is replaced by a limiting density distribution [dV /V) / dr as shown in Fig. 2. The 
density distribution's discontinuities reflect a preferred set of short closed geodesies 
and may be the subject of a future paper, but nevertheless hold no importance for 
the present study. 



''One may substitute the last scattering surface at 2 ~ 1100 for the absolute horizon at 2 = oo 
with little effect on the horizon radius. 

''Given a point p in a multiconnected space, the injectivity radius at p is defined to be the radius 
of the largest ball centered at p whose interior does not overlap itself. Equivalently, twice the 
injectivity radius is the length of the shortest topologically nontrivial path from p to itself. 
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Fig. 1. An injectivity profile shows how much of a manifold's volume has a given injectivity 
radius. 
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Fig. 2. As the bin width goes to zero, the histogram becomes a density distribution. 



2.2. Computational Overview 

To construct an injectivity profile, begin with a fundamental domain for the mani- 
fold. The present study used a Dirichlet domain computed by the computer program 
SnapPea.u SnapPea presents the Dirichlet domain in the Klein model of hyperbolic 
space, where it looks just like an ordinary Euclidean polyhedron. Cover the Klein 
model with a grid of points. The grid is rectangular relative to the Euclidean ge- 
ometry of the model, although of course it's nonrectangular relative to the intrinsic 
hyperbolic geometry. The present study used a 200 x 200 x 200 grid. The program 
scanned the grid, first rejecting grid points lying on or beyond the sphere-at-infinity 
in the Klein model, and then rejecting points lying outside the Dirichlet domain. 
For each grid point lying within the Dirichlet domain, the program computed the 
injectivity radius at that point (Subsection |2.3| ) to determine the correct bin in the 
histogram, then computed the true hyperbolic volume of the surrounding grid cell 
(Subsection 2.4) and added that volume to the correct bin. After all grid points 
were processed, the program printed the volume in each bin to a file. 



2.3. Computing the Injectivity Radius 

At first glance, computing the injectivity radius at a point p is trivially easy: just 
apply elements 7 of the holonomy group P to p and see what the minimum trans- 



2, 2008 1:46 WSPC/Guidelines-IJMPA NearlyFlatHyperbolic 



4 Weeks 

lation distance dist(p,"f{p)) is^ The injectivity radius will be half that minimum: 
^injip) = min^gr ^^^^^^^^^IeII ^ The only problem is that T is an infinite group, so 
we must decide ahead of time how many elements — and which elements — of F to 
apply. 

Our selection criterion for choosing a finite subset of F depends on the con- 
cept of an isometry 7's basepoint translation distance, defined to be the distance 
dist{0,j{0)) that 7 translates the origin O. Say we want to find all 7 G F that 
translate a given point p a distance at most £. Any such isometry 7 has basepoint 
translation distance dist{0,j{0)) < dist{0,p) + dist{p,j{p)) + dist{-f{p),j{0)) < 
Dout + £ + Dout, where Dout is the outradius of the Dirichlet domain D. In other 
words, if we consider all isometries 7 with basepoint translation distance at most 
£ + 2D out, we are sure to find all translates oi p closer than the distance t. 

How large must £ be? A simple argument^ shows that every point p has a 
translate ^{p) lying at a distance less than 2_Di„ + 2D out- So if we were to choose £ 
in the preceding paragraph to be 2D in + 2D out , we would be guaranteed to find the 
nearest translate of every point p, and would therefore know the injectivity radius 
Tinjip)- This method, while rigorous, is inefficient: the number of isometries 7 with 
basepoint translation distance at most £ grows exponentially in £, and the majority 
of those isometrics are unneeded. A more efficient algorithm is simply to guess a 
plausible value of £, find all isometries 7 with basepoint translation distance at most 
£ + 2Dout, and verify afterwards that for each point p we found a translate 7(p) 
at a distance dist{p,j{p)) less than £. There is no guarantee that such a translate 
will always be found, but if one is found then we rigorously know the minimum 
value of dist(p,"/{p)) and hence the injectivity radius rinj{p) — dist{p,^(p))/2. In 
the present study a value of ^ = 1.3 worked in all cases. 

Given a Dirichlet domain D and its face pairings, how do we find all isometries 
7 G F with basepoint translation distance dist{0, "f{0)) at most £7 A simple recur- 
sion does the job: start with the the original Dirichlet domain D; add its immediate 
neighbors in the universal covering space, which we think of as the translates j{D) 
for each of D's face pairing isometries 7; then add the neighbors' neighbors, which 
correspond to products 7271 (-D) of pairs of face pairing isometries; and so on. Con- 
tinue recursively, keeping those images 'y{D) wit]ijlist{0,j{0)) < £ and discarding 
those with dist{0,^{0)) > £. Proposition 3.1 of shows that this algorithm finds 
all isometries with dist{0,j{0)) < £; that is, we needn't worry that the recursion 
will terminate when it encounters only unwanted translates with dist{0, 7(0)) > £, 



■^The holonomy group F is also known as the group of covering transformations. We follow the 
convention of saying holonomy group when thinking of F geometrically as a group of isometries, 
but saying group of covering transformations when thinking of F topologically as a group of 
homeomorphisms. 

■^Let 7 be the face pairing isometry taking one of the Dirichlet domain's faces closest to the origin 
O to its mate. Because the chosen face is maximally close to the origin, distiO ,^(0)) = 2Z)i„, 
where Din is the inradius of the Dirichlet domain D. Thus for any point p in the Dirichlet domain, 

dist{p, 7{p)) < dist{p, O) + dist{0, 7(0)) -I- dist{-y(0), 7(p)) < Dout + 2Di„ + Dout- 
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■ 3. The Klein disk projects radially onto the hyperboloid model of hyperbolic space, taking 
a sphere of Euclidean radius r in the Klein disk to a sphere of intrinsic hyperbolic radius p. 



while some wanted translate with dist{0,"f{0)) < i remains hidden beyond them. 



2.4. Computing the Volume of a Grid Cell 

Inscribe the Klein model of hyperbolic space in a cube with corners at (±1, ±1, ±1), 
and cover the cube with a 200 x 200 x 200 grid. Each grid cell is then a small cube 
of side length 10~^ and Euclidean volume 10~®. The hyperbolic volume represented 
by each grid cell differs from the apparent Euclidean volume, and varies from cell 
to cell. 

It's easy to compute the ratio of the true hyperbolic volume to the apparent 
Euclidean volume. The grid cells are small enough that the volume ratio doesn't 
vary much within a given cell, so pick any point (a;, y, z) within the cell, and let 
r = \J -\- ip' -V be its distance from the origin. By symmetry the volume ratio 
depends only on r, not on the particular point (x, y, z). So take a spherical shell of 
radius r and thickness dr, and compare its true hyperbolic volume to its apparent 
Euclidean volume. The apparent Euclidean volume is 

Ve = 47rr2dr. (1) 

To compute the shell's true hyperbolic volume, project from the Klein model onto 
the hyperboloid model of hyperbolic space (Fig. 3). Similar triangles show that 
the hyperbolic radius p has coshp = \j\J\ — and sinhp = r /\/\ — . Thus the 
area of the sphere is 47rsinh^/o — 47rr^/(l — r^). To find its thickness, note that 
tanh p = r so dp = cosh^pdr = dr / {1 — r^). Hence the hyperbolic volume of the 
spherical shell is 

Vh = 7—-^ (2) 



and the required ratio is 

Vh 4:nr^dr/{l 



(3) 
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3. Horizon Radius 

The variable rhor represents the horizon radius in comoving coordinates, that is, 
in units of the curvature radius of the ambient hyperbohc space. In a standard 
Friedmann-Lemaitre model with matter density today Qm and cosmological con- 
stant f2A, one may compute rhor by integrating the path of a photon backwards 
from z = to z = cxD, obtaining 

dz 



rhor = \/\^m + ^^A - 1| / , (4) 

^ Jo + HaKz + 1)2 + + 1)3 

Replacing the upper limit of z = oo with the redshift z = 1100 of last scatter would 
have little effect on the result. 

For comparison with Gomero et al.p the remainder of the present article will 
consider two sets of plausible choices for the density parameters: 

• r2,„ ~ 0.37 and Q\ = 0.61, for which ^totai = 0.98 and rhor — 0.43, and 

• Urn — 0.37 and JIa — 0.62, for which iltotai = 0.99 and rhor — 0.30. 



4. Results 

A computer program, written in the C programming language and freely available 



from ftp : //ftp .northnet . org/weeks/NearlyFlatHyperbolic, computes injec- 



tivity profiles using the algorithm outlined in Subsection 2.2 along with the details 
explained in Subsections ^.3| and |2.4 For the first ten low- volume hyperbolic mani- 
folds from the Hodgson- Weeks census,til which include the seven manifolds consid- 
ered by Gomero et aL,u the program obtains the injectivity profiles shown in Fig. 4. 
Shaded backgrounds mark the cutoffs rhor = 0.43 and rhor = 0.30, corresponding 
to ^Itotai = 0.98 and ^Itotai = 0.99, respectively. Only in the fraction of the man- 
ifold to left of each cutoff would the topology be detectable. Table 1 summarizes 
the results. 



Table 1. The fraction of each manifold's vol- 
ume in which its topology is potentially detectable 
when ntotal = 0.98 or 0.99. 

manifold fttotal = 0.98 Utotal = 0.99 
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Fig. 5. Injectivity profiles for the two smallest limiting cusped manifolds. 

The manifolds in Table 1 contain no closed geodesies shorter than 0.36. However, 
there exist families of hyperbolic 3-manifolds containing arbitrarily short closed 
geodesies. The manifolds within each family approach a limiting cusped manifold 
which is, in effect, a manifold with a geodesic of length zero. Furthermore, the in- 
jectivity profiles of the manifolds in each family approach the injectivity profile of 
the limiting cusped manifold. Fig. 5 shows the injectivity profiles for the two small- 
est limiting cusped manifolds, mOOS and mff-^,^ while Table 2 shows the fraction 
of each in which an observer could potentially detect the topology. Note that for 
^totai close to one, the fraction of the manifold in which the topology is detectable 
is roughly proportional to 1 — ^Itotai ■ The Dirichlet domains for these cusped man- 



ifolds have infinite outradius, so in the algorithm of Subsection 2.S we considered 
all isometrics of basepoint translation distance less than the plausible but arbitrary 
value of 2.5; thus these results for cusped manifolds are not rigorous, but they are 
most likely correct nonetheless. The beginning of each profile is noisy because at 
most a few grid points fall into each bin; refining the grid diminishes this effect at 
the expense of a slower computation. 

Table 2. The fraction of each cusped manifold's volume in which its topology 
is potentially detectable when i^total is close to one. 



manifold 


^total = 0.98 


^total = 0.99 


^totai = 0.995 ntotai = 0.9975 


m003 


17% 


8% 


4% 2% 


m004 


67% 


32% 


16% 8% 



''The cusped manifolds mOOO, mOOl, and m002 are nonorientable and not the limit of any sequence 
of closed manifolds. 
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5. Conclusions 

In the case iltotai — 0-98, the horizon radius rhor — 0.43 is large enough that 
an observer at a random location in a small hyperbolic universe would have a 
reasonable chance of detecting the topology. However, in the case ^totai = 0.99, the 
horizon radius Vhor drops to 0.30 and a random observer would have little or no 
chance of detecting the topology. 

There exist closed hyperbolic 3-manifolds with arbitrarily short closed geodesies, 
so no matter how close fltotai is to one, there will always be infinitely many hyper- 
bolic manifolds in which well-placed observers could detect the topology. However, 
the probability that an observer could detect the topology from a random point in 
such a manifold decreases in proportion to 1 — ^totai- 
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